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Necessary conditions that the spurious state associated with the translational motion and its 
double-phonon state have zero excitation energy in extended RPA (ERPA) theories which include 
both one-body and two-body amplitudes are investigated using the small amplitude limit of the 
time-dependent density-matrix theory (STDDM). STDDM provides us with a quite general form 
of ERPA as compared with other similar theories in the sense that all components of one-body 
and two-body amplitudes are taken into account. Two conditions are found necessary to guarantee 
the above property of the single and double spurious states: The first is that no truncation in the 
single-particle space should be made. This condition is necessary for the closure relation to be used 
and is common for the single and double spurious states. The second depends on the mode. For 
the single spurious state all components of the one-body amplitudes must be included, and for the 
double spurious state all components of one-body and two-body amplitudes have to be included. 
It is also shown that the Kohn theorem and the continuity equations for transition densities and 
currents hold under the same conditions as the spurious states. ERPA theories formulated using 
the Hartree-Fock ground state have a non-hermiticity problem. A method for formulating ERPA 
with hermiticity is also proposed using the time-dependent density-matrix formalism. 

PACS numbers: 21.60.Jz 



I. INTRODUCTION 

The double phonon states of giant resonances have become the subject of a number of recent experimental and 
theoretical investigations 0,0- I n the case of giant resonances (single phonon states), the random phase approximation 
(RPA) has extensively been used as a standard microscopic theory to study basic properties of giant resonances . It 
is guaranteed in RPA that physical states do not couple to spurious states such as that associated with the translational 
motion because RPA in the Hartree-Fock (HF) basis gives zero excitation energy to spurious states. For a microscopic 
study of the double phonon states of giant resonances, we need to extend RPA to deal with two-body amplitudes as 
well as one-body amplitudes. One of such an extended RPA theory (ERPA) may be the Second RPA (SRPA) .4] which 
has also extensively been used to study decay properties of giant resonances 0, ||| . When the double phonon states are 
studied in ERPA, it should be guaranteed that both spurious states and their double-phonon states are decoupled from 
physical states. The aim of this paper is to investigate necessary conditions that the spurious state associated with the 
translational motion and its double-phonon state have zero excitation energy in ERPA. We use the small amplitude 
limit of the time-dependent density-matrix theory (STDDM) 0. The reason why STDDM is used is that containing 
all components of one-body and two-body amplitudes, STDDM constitutes a more general framework of ERPA than 
SRPA. It will be shown that keeping all components of the one-body and two-body amplitudes in ERPA is essential 
in bringing the spurious states to zero excitation energy. Any ERPA including STDDM , which is formulated using 
an approximate ground state, inherently has asymmetry and non-hermiticity. A method for recovering symmetry and 
hermiticity in the framework of the time-dependent density-matrix formalism is also proposed in this paper. The 
paper is organized as follows: STDDM is presented and its relation to other ERPA formulations is discussed in Sect. 2. 
Necessary conditions to give zero excitation energy to the spurious state associated with the translational motion 
and its double-phonon state are discussed in Sect. 3. The Kohn theorem 0, [tj ^| and the continuity equations for 
transition densities and currents are also discussed as related subjects in Sect. 3. In Sect. 4 a method for formulating 
ERPA with hermiticity is proposed and Sect. 5 is devoted to a summary. 

II. EXTENDED RPA FORMALISM 

A. Small amplitude limit of the time-dependent density-matrix theory 

The time-dependent density-matrix theory (TDDM) gives the time-evolution of a one-body density-matrix p(l, 1') 
and the correlated part C(12, 1'2') of a two-body density-matrix where numbers denote space, spin, and isospin 
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coordinates. Linearizing the equations of motion for p and C with respect to 6p and 5C, where 5p and 5C denote 
deviations from the ground-state values p and Co i.e. dp = p — po and SC — C — Co, respectively, we obtain 
STDDM0. Expanding Sp and 5C with single-particle states ipa as 

Sp(ll',t) = (1) 



5C(12l'2',i) = ]T X a p a ,p,(t)Mht)M2,tWa'(l',tWp(2',t), (2) 

q/3q'/3' 

and assuming the HF ground state, that is, po is the one-body density-matrix in HF approximation and Co = 0, we 
obtain the following equations of STDDM for the Fourier components of x aa >(t) and X a p a ipi(t) Q: 



(<jj - e a + e a /)x aa < = (f a > - f a )y^^(a\\v\a'\') 



AX\>\ 



XX' 



[XxA> a 'X»{a\"\v\\\'}-X aA'AA"(AA"|w|a'A'}], (3) 



AA'A" 



(uj - e a - ep + e a - + ep>)X a p a ,fji = - ^[(fpfa'fp' + fpf a 'fp')(Xf3\v\a'f3') A x aX 

A 

+ (fafa'fp' + f a fa'f/3')(aX\v\a'P') A X/3X ~ (fafpfp' + fuf/3fl3>)(af3\v\X(3') A X Xa > 
- (fatfpfa 1 + faf(jf a '){al3\v\a'\) A Xxfi'] 

+ ]T[(1 - fa - M(a(3\v\XX')X xx , a , , - (1 - f a , - fp,)(XX'\v\a'p')X a p X x] 

AA' 

+ - f a ){aX\v\a'X'} A X x , 0Xf3 , - (f fj , - f a ){aX\v\X' 0') A X x ,p a , x 

AA' 

+ (//3< - fp)(Xp\v\X'P') A X aya , x - (f a , - fp){Xp\v\a'X') A X aX . X p,], (4) 

where e a is the HF single-particle energy, f a — l (0) for occupied (unoccupied) single-particle states and = 1 — /«, 
and the subscript A indicates that the corresponding matrix element is antisymmetrized. Let us mention that eqs.Q 
and (@J may also be obtained from the following equations of motion 

($o|[a>a,ff]|$) = iu(^ \a+a a \^), (5) 
{$o\[a+a+a p a a ,H}\<S>) = u{$> o \a+,a^,a a a \$>), (6) 

where [ ] is the commutation relation, H the total hamiltonian consisting of the kinetic energy term and a two-body 
interaction, |$o) the ground-state wavefunction and $} the wavefunction for an excited state with excitation energy w. 
Linearizing eqs.|0 and I© with respect to x aa ' = (3>o|a^/a Q |<l>} and X a p a ip> = ($>Q\a~^,ap,apa a \$>), and assuming the 
HF ground state for |$ ) when expectation values for the ground state are evaluated such as ($o|ai a Q' |$o) ~ 
we obtain eqs.® and (@J). 

In the following we discuss some relations of STDDM with RPA and other versions of ERPA. When the coupling 
to the two-body amplitudes X a fj a ipi is neglected in eq.®, the equation for the one-body amplitudes becomes 

{uj - e a + e a ,)x aa > = (f a > - f a )^2(aX\v\a'X') A x x > x . (7) 

AA' 

When f a is the Fermi-Dirac distribution, eq. is equivalent to the finite temperature RPA 0> EH • Hereafter single- 
particle indices p and h are used to refer to unoccupied and occupied single-particle states, respectively. Since the 
sums on the right-hand sides of the equations for x p h and X} lp are unrestricted, x p h and Xh P can couple to x pp i and 
Xhh' ■ Such a coupling scheme of eq. (UJ) may be better understood in matrix form 



€ p -c h + {ph'\v\hp') A (pp'\v\hh') A (pp'\v\hp") A {ph'\v\hh") A \ 

-(hh'\v\pp') A e h -e p - (hp'\v\ph') A ~{hp'\v\pp") A ~(hh'\v\ph") A \ 

e p - t p i I 

e h -e h , J 



/ x p > h > \ 




/ Xph \ 


Xh'pi 


-„ 


Xhp 


1 Xp'ipi 




x pp > 1 

V X h h' ) 


\ Xh"h' ) 





(8) 
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where obvious summation symbols and Kronecker's (Ts are omitted for simplicity. Since the hamiltonian matrix is 
non-hermitian, x aa > is orthogonal not to x aa i but to a left-hand-side eigenvector x aa i which satisfies 

(uj - e a + e a ,)x* aa , = ^2(f\> - f x ){\a'\v\\'a) A x* xx , 

AA' 

= ^^((pa' \v\ha) A x* ph — (ha' \v\pa) A x* hp ) (at temperature T = 0). (9) 

ph 

The matrix form of eq.© becomes 



/ ~* ~* ~* ~* \ 

I Xp'i h" t Xh"p" j Xppi , Xfofoi ) 



( ep - e h + (p"h\v\h"p) A (p"p\v\h"h) A (p"p\v\h"p') A (p"h\v\h"h') A 

(h"h\v\p"p) A e h -e p - (h"p\v\p"h) A -(h"p\v\p"p') A -{h"h\v\p"h') A 
e P -e P ' 



V o o o 

= w \%ph> x hp> x pp'>®hh') • (10) 
The ortho-normal condition is written as 

(A|A') = ^^ a ,(A) w (A') = ^v, (11) 

aa' 

where |A) represents an eigenvector x aa i with the eigenvalue lu x , and |A) the left-hand eigenvector of the hamiltonian 
matrix with the eigenvalue uj x . The completeness relation becomes 

^|A)(A| = ^ W (A)^(A)=7, (12) 

A A 

where I is the unit matrix. These ortho-normal and completeness relations are generalizations of the RPA ones. Due 
to the occupation factor f a — f a i the one-body amplitudes x pp : and Xhw vanish unless uj = e a — e a r (see eq.lJTJ)) 
whereas x aa > always have all components as seen from eq.©: x aa i corresponds to the generalized RPA amplitude 
which appears in the Landau's expression for the damping width of zero sound [bH. If the particle (p) - particle 
(p) and hole (h) - hole (h) components of glected, eq.Q is reduced to the RPA equations, 

(oj - e p + e h )x ph = y\(ph'\v\hp') A Xp> h ' + (pp'\v\hti) A x h >p>], (13) 

p' h' 

(oj - e h + e p )x hp =- }J(hh'\v\pp r ) A Xp> h > + (hp'\v \ph') A x h > p >]. (14) 

p'h' 

When the coupling to the one-body amplitudes is neglected in eq.(@J, the equation for the two-body amplitudes 
become 



(w - e a - ep + e a , + e^X^p, = ^[(1 - f a - fp)(ap\v\\X')X xx , a , p , - (1 - f a , - fp>)(\\'\v\a' p')X a p XX ,} 

AA' 

+ 5Zl(/«' - fa){a\\v\a'\'} A X x ,p k p, - (fp - f a ){a\\v\\'p') A X x ,p a , x 

AA' 

+ (//?' - f )(\0\v\\'P'} A X aX , a , x - (f a , - f ){X/3\v\a'X'} A X aX , X0 >}. (15) 

This equation is equivalent to the formula given in Ref.0] for the two-body space. Keeping only the 2p-2h, 2 h-2p 
and lplh-lplh components of X a p a ip* in ea . leads to the version of ERPA for low-lying two-phonon states [l7j . 
It has been pointed out |17| that the lplh-lplh components of X a p a ipi are important to reproduce collectivity of 
low- lying double-phonon states. A time-dependent version of ea. (|15fl has been applied to the double-phonon states 
of giant dipole and quadrupole resonances in 40 Ca using a realistic Skyrme-type interaction for both the mean-field 
potential and the residual interaction, and it was found that the 2p-2h, 2h-2p and lplh-lplh components are the 
most important two-body amplitudes for these double-phonon states 0. In eqs.© and l@J the one-body amplitude 
x aa' and the two-body amplitude X a p a i p have all components: For example, x aa > has lp-lh, lh-lp, lp-lp and lh-lh 
components. On the other hand only the lp-lh and lh-lp components of x aa i and the 2p-2h and 2h-2p components 
of X a p a ipi are taken into account in SRPA Q and the SRPA equations are obtained from eqs.© and (0J by keeping 
only these amplitudes. 
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Equations JJjJ and Q have asymmetric couplings between the x aa > and X a p a ipi amplitudes: In eq.J2J) x aa ' couples 
to all components of X a p a >p>, while in cq.Q only the 2p-2h, lp-3h and lh-3p components of X a0a ipi (and their 
complex conjugates) couple to x aa i due to the occupation factors (fpfa'fp' etc.). Equations and lfT5)l which 
have no coupling between one-body and two-body amplitudes are also non-hermitian due to occupation factors such 
as f a ' — f a - The asymmetry and non-hcrmiticity originate from the structure of the equations for the reduced 
density matrices (see eqs.(0 and ©)• For a non-hermitian hamiltonian matrix, the left-hand-side eigenvectors of 
the hamiltonian matrix constitute a basis which is orthogonal to (x aa i, X a p a ipi), and the ortho- normal condition is 
written as 

<A|A'} = $>^(A):w(A')+ Kp a ,p,(\)X a0a , l3 ,(\') = S xx ,, (16) 

aa' a(3a' j3' 

where |A) represents an eigenvector (x aa t, X a p a ipt) with the eigenvalue uj\, and |A) the left-hand-side eigenvector of 
the hamiltonian matrix with the same eigenvalue. The completeness relation is written as 

The asymmetry and non-hermiticity in eqs.© and Q are necessary to prove the properties of the spurious states 
and the Kohn theorems as will be discussed below. However, due to the non-hermiticity of the problem some of the 
eigenvalues may become complex. Our exploratory numerical calculations for the oxygen isotopes 22 ' 24 using the 
neutron 2s and Id states and a pairing-type residual interaction which had been used in the calculations of quadrupole 
states in these nuclei [l9j show that the non-hcrmiticity of STDDM is quite moderate: Only a small fraction (about 
10%) of the eigenstates have complex energies, whose imaginary parts are less than 0.1 MeV. The results of these 
numerical calculations will be discussed elsewhere. On the other hand we will show in Sect. 4 that there is a 
prescription for constructing ERPA with symmetry and hermiticity using a correlated ground state in TDDM. 

III. SINGLE AND DOUBLE SPURIOUS STATES 
A. One-body and two-body operators for the translational motion 

We consider the following one-body and two-body operators associated with the translational motion: 

P=^2(a\-iV\[3)a+ap (18) 

a/3 

and 

P-P = J2(a\ - V 2 \a')a+a a , + ^ (a\ - iV\a'} ■ {/3\ - iV\/3')a+ap~a >a a > . (19) 

aa' a{3a' f3' 

Since the hamiltonian H has translational invariance, these operators commute with H, that is, [H, P] = [H, P-P] = 0. 
We will evaluate w($o|P|3?i) and u>{$>o\P-P\$2), where |$i) and |$2) are the spurious states excited with P and PP, 
respectively, and show that these states have zero excitation energy in STDDM. In the case of the exact problem, it 
is, with eqs.® and ©, trivial to see that uj(Qq\P\$i) and lo(<&q\P ■ -Pl"!^) arc identical to zero because the left-hand 
sides of eqs.(|5l and © are reduced to the commutation relations between the hamiltonian and these translational 
operators. Since the linearization and the HF assumption are made in the derivation of STDDM, it is not trivial to 
show the above properties of the spurious states. However, the linearization should be valid in the weak coupling limit 
and therefore we can anticipate that the Goldstone theorem also holds in this case, provided that the linearization 
procedure is correctly performed. 

B. Spurious states in RPA 

RPA gives zero excitation energy to the spurious state |$i) excited with P, although only the lp-lh and lh-lp 
components of the one-body amplitudes are taken into account. To illustrate our approach for the problem of the 
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spurious states, we begin with proving w($o|-P|^i) = in RPA. Using the relation (^ol^a'^l^i) — x aa< and eas. H13() 
and IT4*jl for x p h and Xh P , we modify uj($q\P\$i) as 

w($o\iP\$i) = u^2[(h\V\p)xp h + (p\V\h)x hp \ 

ph 

= ^2[(h\V\p)(e p - e h )x ph + (p\V\h)(e h - e p )x hp ] 

ph 

+ [(h\V\p)({pti\v\hp')AX P 'h' + (pp'\v\hh') A x h > pl ) 

php'h' 

~ (p\V\h)((hp\v\pti) A x h > p > + {hti\v\pp') A )x P >h']- (20) 

A further modification is made using hoip a — e a tp a , where ho is the HF single-particle hamiltonian, and the closure 
relation ip p (r)ip*(r') = S 3 (r- r') - Y,h^h{r)^l{r'): 

w{$ |iP|$i) = ^[{h\\V ,ho]\p)x P h + (h\[V,h ]\h)x hp \ 

ph 

+ ^[{{hti\V lV \hp) A - (hh'lvV^hp) + (hh'\vV 2 \ph))x ph/ 

hph' 

+ {{hp\Viv\hti} A -{hp\vVi\hti) + {hp\vV 2 \tih})x h , p )}. (21) 
The first term on the right-hand side of the above equation can be written in terms of v using 

(a'\[V,ho]\a) = ^[(a'ft|(Viw)|a/i)^ 

h 

- {a'h\vVi\ha) + {a'h\vV 2 \ha)}. (22) 

Finally ea. 121(1 becomes 

w($oKP|$i) = J2[({hh'\{Viv)\ph') A - (hh'\vVi\h'p) + {hti \vV 2 \ti p))x ph 

phh' 

+ ((ph'\(Viv)\hh') A - (ph'\vVi\tih) + {ph'\vV 2 \h'h))x hp ] 
+ ^[{{h'h\{Vxv)\h'p) A - (h'h\vVi\pti) + {h'h\vV 2 \ph'))x ph 

phh' 

+ {{h'p\(V 1 v)\h'h} A -{h'p\vV 1 \hh'} + {h , p\vV 2 \hh , ))x hp )}, (23) 

where (Vi«) means that Vi acts only on v. Since v has translational invariance, the sum of the following two terms 
on the right-hand side of eci. ((23(l becomes zero 

(hh'\(Vxv)\ph') A + (h'h\(Viv)\h'p) A = (hh'\(y lV ) + (V 2 v)\ph') A = 0. (24) 

Another sum of the two terms also vanishes 

- {hh'\vVx\tip) + {h'h\vV 2 \ph') = (h'h\ -vV 2 +vV 2 \pti) = 0. (25) 

Similarly, all other terms on the right-hand side of ea. 1(23(1 cancel out. This means co = 0. As shown above, both the 
inclusion of the backward amplitude Xh p and the unrestricted sum over unoccupied single-particle states are essential 
in RPA to give zero excitation energy to the spurious state. 



C. Spurious state in STDDM 

Along the lines illustrated above, we then show that w($o|-P|$i) = in STDDM. Using the equation for x aa i 
(eq.©), we modify uj($q\P\$i) as 

w($ |iP|$i) = ^2{a'\V\a)cJx aa , 

a a' 
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= y^(a / |V|a){(e a - e a >)x aa < 

olol' 

+ (fa' ~ fa) y^^aX\v\a'\') A xy X 

XX' 

+ X) [Xxx>a>x'>(aX"\v\XX')-X a x'xx»(XX"\v\a'X')}}, (26) 

AA'A" 

where the sums are over both occupied and unoccupied single-particle states. Using hoip a — e a ip a and the closure 
relation ip(r)ip^(r') = 5 3 (r — r'), we further modify eq. (|26|l as 

w($o|*^|*i) = Y,(a/\\y,ho]\<*)**a> 

aa' 

+ [(hX\Viv\h\')A - ((h\\vVx\hX') - (hX\vV 2 \X'h))}x X 'x 

XX' h 

+ [Xxx'a'X"(a'X"\\7 1 v\XX')-X a y XX „(XX"\v\7 1 \aX')}. (27) 

aa' XX' X" 

Using ea. H22(l and changing summation indices, we obtain 

w($ |iP|$i) = ^2[(a'h\(Vxv)\ah}A-{a'h\v\/i\ha) + {a'h\vV 2 \ha) 

aa' h 

+ (ha \(X7 iv)\ha) a — (ha'\v\7\\ah) + (ha \vV2\ah))]x aa i 

+ ]T X*/3a'0'(a'0'\(Viv)\aP). (28) 

a0a'f}' 

The first sum on the right-hand side of the above equation which includes terms proportional to x aa i is a generalization 
of ea. l|23|l and vanishes for an interaction v with translational invariance. The second term on the right-hand side of 
eq.J2HJ can be expressed as 

X a(3al0 ,(a'p'\(S7iv)\ap) 

afia'P' 

= ~ X a p a ,p,({a'fi'\(V 1 v)\a[}) + (l3'a'\{V x v)\l3a)) 

a/3a'(3' 

= \ X a0a>pi (a'P'\(V 1 v) + (V 2 v)\a(3). (29) 

a/3a'/3' 

Since v has translational invariance, (Viv) + (\7 2 v) = 0. Thus again uj(&o\P\&i) = is proven. As shown above, 
unrestricted summation over single-particle indices a and a' in ea. il26|) is essential to derive the last term on the 
right-hand side of eq. (|27|l . This means that any ERPA formalisms with restricted one-body amplitudes cannot give 
zero excitation energy to the spurious state associated with the translational motion. However, this does not depend 
on approximations for two-body amplitudes as long as the symmetry property is respected as seen in eq. <|29|) . 

D. Double Spurious state in STDDM 

In a way similar to the above, we show that uj($o\P ■ P| ( E > 2) = 0, where |<I>2) is the double spurious state. The term 
w($o|-P • P\® 2 ) contains both the one-body and two-body contributions, 

- lu(^ \P ■ P\<S> 2 ) = ^{^((a'|V 2 |a) -^2<a'|V|/i) • (h\V\a))x aa , 

aa' h 

+ ( a 'M a ) ■ <P'W0)Xalla>P>}. (30) 

a(3a'/3' 
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Using eqs.© and Q for x aa > and X a p a ipi, h ip a = £ a ipa and the closure relation J2 a ?/>(r*)?/>* (r') = 8 3 (r — r'), we 
modify the right-hand side of the above equation. After some lengthy manipulations, the terms containing x aa i and 
one summation index over occupied single-particle states become 

2 [{ot'h\{V\v)\ah) A + (a'h\(Vxv) ■ Vi\ah) 

aa'h 

- (ah\(Viv) ■ Wi\ha) + {ha\(Viv) ■ Vi\ha) - (ha'\(\7 lV ) ■ Vi\ah)]x aa , 
+ 2 £ [(a'h\(Vi ■ V 2 v)\ah) A + (a'/i|(Viv) ■ V 2 \ah) 

aa'h 

- {a'h\{V lV ) ■ y 2 \ha) + (a'h\(V 2 v) • Vi\ah) - (a'h\(y 2 v) • Vi\ha}]x aa ,, (31) 

where the first sum comes from the terms with x aa > on the right-hand side of ea. l|30(l and the second sum from 
the term with X a p a ifji. Since v has translational invariance, (S7 2 v) + (Vi • ~S/ 2 v) = 0. Therefore, the sum of the 
following two terms in eq.JjHJ becomes (a' h\(V 2 v)\ah) a + (a'h\(Vi ■ V 2 v)\aK)A = 0. All other terms vanish for 
similar reasons. In addition to the terms shown in ea. (|31|l . there appear terms with x aa > and two summation 
indices over occupied single-particle states, and also terms with X a {j a ipi in the modification process of ea. H3()(l . It is 
straightforward, though lengthy, to show that these terms also vanish for a translationally invariant interaction. Thus 
uj($>o\P • P\®2) = 0, that is, uj = 0. As mentioned above unrestricted summation over single-particle states is again 
essential to obtain this conclusion. This means that only ERPA's with all one-body and two-body amplitudes, that 
is, Xpii , x/ip, Xpp' , Xfiii' , Xpp'hfi' , Xhh'pp' , Xphp'h' , Xphwh" , Xh>fi"ph , Xpp'hp" , Xhp"pp' , Xpp'p"p>» find ^/i'/t"/i'"i give zero 
excitation energy to the double-phonon state corresponding to the spurious mode associated with the translational 
motion. 



E. Single and double Kohn modes 

When a system is confined to a harmonic potential U = Imw^r 2 , the spurious mode associated with the translational 
motion has an eigenvalue of hujQ, inde pen dently of the translationally invariant two-body interaction. This property 
is known as the Kohn theorem 0, 0, llf)j | . In this subsection we show that our ERPA equations satisfy the Kohn 
theorem and also that the eigenvalue of the double Kohn mode becomes 2Tiuiq. Due to the presence of the harmonic 
potential the single-particle states are chosen to be eigenstates of the modified hamiltonian, h'ip a = e a 4>a, where 
h' = ho + ^mujgr 2 . In a way similar to the spurious mode, we evaluate w(<f>o|-P|^ ) i} using the equations of motion in 
ERPA. Since the two-body interaction has translational invariance, terms with the two-body interaction vanish and 

w($o|*£|*i) = - ^{oi\mwlr\a')x a , a = -mwg($ |3|*i) (32) 

act' 

holds, where Q = ^{a^a^a^aa' . Similary, non- vanishing contribution to o;( < I>o|Q| ( i ) i) comes only from the kinetic 
energy term, and we obtain 

h 2 

w($ |Q|$i) = ($ |iP|*i>- (33) 

771 

It is essential to keep all components of the one-body amplitudes to obtain the above expressions. From eqs.J23) an d 
(|33fl . we get lj = ±?iu;o- 

In the case of the double Kohn mode, expectation values of three operators couple in the following way, 

u;(^P-iP\^2) = 2mw 2 ($ \Q-iP\$ 2 ), (34) 

h 2 

w($ Q |Q-£P|$ 2 ) = — ($ Q \iP-iP\$ 2 ) + mujl($ \Q-Q\<f> 2 ), (35) 

777 

- - h 2 - - 

W ($ |Q-Q|$ 2 ) = 2— ($ |Q.iP|*a). (36) 

TO 

The right-hand side of ea. H34l) comes from the harmonic potential. Both the kinetic energy term and the harmonic 
potential contribute to the right-hand side of eq. I|35() , and the kinetic energy term becomes non- vanishing on the right- 
hand side of ea. l|3t)|) . All terms with the two-body interaction vanish due to translational invariance. It is essential 
to keep all components of the one-body and two-body amplitudes to derive eas. (|34ll - (|36|l as in the case of the double 
spurious mode discussed in Subsect. 3.4. From the above equations we get u> = ±2Tuvo- 
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F. Continuity equations 

We end this section by showing that our ERPA equations satisfy continuity equations. In a way similar to the 
single spurious mode we evaluate uj\($>o\p(r)\&\) , where p is the density operator p(r) = ^ ip^,(r)tp a ' (r)a^a a i , and 
obtain 

wa<*o|#OI*A> = -V • ($o|j(r)|$A), (37) 
where the current operator j is given by 

^ = ^ Z>a(*DWa,(f) - (VC(r1)^(r1]a+a a , (38) 

for a momentum-independent two-body interaction. Thus the continuity equation for the one-body transition density 
and current is satisfied. Keeping all components of the one-body amplitudes is essential to obtain the continuity 
equation. 

Similarly, the transition amplitude for the two-body density operator p^if^r 1 ) defined by 

p 2 {f,r)= ^a{^i>p(^)^P'(^)i>a'(^ a i a p a P' a c'' (39) 

satisfies the continuity equation 

Ux($o\P2(?,7)\$x) = -(V r • ($o|J 2 (r,r')|$ A ) + V J . • ($ |j 2 (^, »DI*A», (40) 
where the two-body current operator j% for a momentum-independent two-body interaction is given by 

J 2 (ry) = |;E[CM(Vi'«) - (VCtrDi-frDl^t^t^a+a+a^^. (41) 
In the derivation of ea. l4U|) it is again essential to keep all components of the one-body and two-body amplitudes. 

IV. ERPA WITH HERMITICITY 

The equations of STDDM (eqs.lJSJ and Q) show asymmetry and non-hermiticity, although this causes no problem 
in conserving various physical properties as discussed above. In the following we show that ERPA with symmetry 
and hermiticity can be formulated using the equation-of-motion approach |20j | and the correlated ground state in 
TDDM. We have pointed out (2l| , in deriving the Landau's expression for the spreading width of a collective state, 
that it is important to include ground-state correlations to remove the asymmetry in STDDM. It is well-known [2(| 
that the asymmetry problem always appears in the equation-of-motion approach when the ground state is replaced 
by an approximate one. Before presenting the formulation of our ERPA, therefore, we summarize the origin of the 
asymmetry in the equation-of-motion approach. When |$o) is the exact ground state of the hamiltonian, there exists 
an identity involving a one-body operator A = a+a a i and a two-body operator B — a~^a~^ap>a a i: 

($ \[[B,H],A]\$o) - ($o\[[A,H],B]\$o) = (®o\[H, [A,B\]\* ) = 0. (42) 
When |$o) is approximated by the HF ground state, the above identity is violated, that is, 

($ \[H,[A,B]]\*o) t^O (43) 

and, consequently, 

($o\[[B,H],A]\$ Q ) + ($o|[A,ff],B]|<S ). (44) 

Since the left-hand side of the above equation describes the coupling of the one-body amplitudes to the two-body ones, 
and the right-hand side, that of the two-body amplitudes to the one-body ones, the resulting ERPA has asymmetric 
couplings. In order to avoid the difficulty of ea . (|44|l . Rowe introduced a symmetrized double commutator |20| . 
However, it was pointed out jl5| that there is an ambiguity in the choice of such a double commutator. 
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Now we proceed to the presentation of our ERPA with ground-state correlations. The ground state |$o) in TDDM 
is constructed so that 

($o\[H,a+a a ,]\$ ) -0 (45) 

and 

($o\[H,a+a^af3,a a ,}\$ ) = (46) 

are satisfied for any single-particle indices 7] . In other words the occupation matrix n° a , and the correlation matrix 
Cafia'p'' * ne expansion coefficients of po and Co, respectively, are determined in TDDM so that the above two equations 
are satisfied. The explicit expression for eas. (|45|l and l|46|l depends on the single-particle state ip a . The equations 
for -nP aa , and C°,^, Q/3 shown in Appendix A are obtained when ip a is chosen to be an eigenstate of the mean field 
hamiltonian ho(po), that is, 

n 2 v 2 r 

ho{po)i) a {l) = — ^-^W + / d2«(l,2)|>o(2,2)^ a (l)-po(l,2)V»a(2)]=eaV'a(l), (47) 

where 

po(H') = E<^U)C'(1')- (48) 

act' 

Although it is not evident to find an analytic solution of eas.l|45|l and l|46(l |2^ . a method for obtaining n° aa , and 
Capa'p' numerically has been proposed [23l | and already been tested for realistic nuclei in the study of giant resonances 
built on the correlated ground state 0,E1]- Since the commutation relation [A,B] = [a+a a i, a^a^aya^] in ea. ll4*2|) 
becomes a sum of two-body operators, we find 

($o\[H, [a+a a/ ,a+a+a Y ap,}}\<S>o) =0 (49) 

which holds due to ea. l|46[l . This means that the coupling matrices are symmetric, that is, 

{^^[[a+a+a^ap, , ff], a+a a >]|$ ) = (^o\[[a^a a/ ,H],a^a+a Y a^]\^ ) (50) 

for the correlated ground state in TDDM. The ERPA equations based on the TDDM ground state are formulated 
using the equation of motion approach [2£j as 

{y Q \{[a+a a ,,H},Q+}\y Q ) = oj{y Q \[a+a a ,,Q + }\y Q ) (51) 

(y Q \[[:a+a+ap,a a , :, H], Q + ]|*o) = w<*o|[: a+a+a^a^ :,Q + ]|* ), (52) 

where the operator Q + is defined by 

Q + = ^2(x\ya+a X ' + ^ \ 1 \ 2 x' 1 \' 2 ■ a^a^a^a^ :) (53) 

and |^o) is assumed to have the following properties 

Q+\*o) = m (54) 
Ol^o) = 0. (55) 

In eqs.JHSJ an d ((HSJ, : : stands for : a+a~papia a i := a^a^a^aa' — A(a+a a > {^^a^api \$>o) + <ipapi {^o\a+ a a > |^o)), 
where A is an antisymmetrization operator. The above equation can be written in matrix form 

A C \ ( x \ ( S\ T\ \ I x 



b Djyxj'^yr, s 2 )\x)> ^ 

where each matrix element is given by 

Stia'a-.XX') = (* |[a+a Q ,,a+ aA ']|*o}, (57) 

S 2 (a'f3'af3 : AiA 2 A' 1 A^) = (tf„|[: a+a+a^a^ :, : a+a+ 6^0^ :]|* ), (58) 
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T^a'a : A 1 A 2 A' 1 A' 2 ) = (V \[a+a a ,, : a+ a+ a A >a A; :]|* ), (59) 

T 2 (a'P'a(3 : AA') = (* |[: a+a+a^ :, 4°v]|*o>, (60) 

A(a'a:AA') = (tf |[[a+<V, fl], «a°A']|*o}, (61) 

B(a'f3'af3 : XX') = (tf„|[[: a+a^/a^ :, ojovll^o), (62) 

C(a'a : A 1 A 2 A' 1 A' 2 ) = (tf |[[a+a a ,, fl], : « ^ a A / :]|*„>, (63) 

D(a'f3'af3 : XM^) = (*o|[[: a+a+a^a,,, :, : a+a+ a ^oaj :]|*o>. (64) 

When the above matrices are evaluated, the ground state j^o) is replaced by |<&o) m TDDM. Then all matrices in 
the above are written in terms of rfi aa , and 0^g a ,g,, which are shown in Appendix B. Due to eas. 1)45(1 and ((46(1 . the 
above matrices have the following symmetries 

A(a'a : XX') = A(X' X : aa') = A(AA' : a'a)*, (65) 
B{a'f3'af3 : AA') = C*(A'A : a(3a'/3') = C(XX' : a'/3'af3)*. (66) 

This version of ERPA gives zero excitation energy to spurious modes associated with operators O which commute 
with H and consist of one-body and (or) two-body operators. This is because w(5'o|C| x I') = (T'o | [-ff, C] = 
holds due to eas. (|51|l and Although the coupling matrix between the one-body and two-body amplitudes is 

symmetric, the hamiltonian matrix on the left-hand side of eo. 1(56(1 is not yet hermitian because D = D + does not 
hold. This originates in the fact that ($o|[-ffj [■ o,^ataf3'a a ' :, : a^a^ay^ax^ :]]|$o) 0- In order to obtain a hermitian 
hamiltonian matrix without any truncation of the two-body amplitudes, we need to impose 

($ |[#,ai4 a 7 a 7' a /3' a "']l $ o) = (67) 

in addition to eas. 1(45(1 and ((46|) . This condition guarantees {$q\[H, [: a^atapra a / :, : a^a^ctx' a\> :]]|$o) — 0, and 
thereby 

D(a'(3'a(3 : X 1 X 2 X' 1 X' 2 ) = D(A' 1 A 2 A 1 A 2 : a/3a'/3') = D(X 1 X 2 X' 1 X' 2 : a'f3'a(3)*. (68) 

Equation 1167(1 is explicitly shown in Appendix B. For a hermitian hamiltonian matrix the ortho-normal condition is 
given by 26 

K'K'){% £)(£;) =w, (69) 

where x^ and A M constitute an eigenstate of ea. ((56(l with u) = ui^. The completeness relation becomes 

The transition amplitudes for one-body and two-body operators, z = (^o\a^a a i | M/) and Z = (*&o\ '■ <J'a a p a f3' a a' '■ 
respectively, are calculated as follows 



z \ _ ( Si Ti \ ( x 

z -\t 2 s 2 \ x 



(71) 



Equation l(56|) has a certain similarity with the so-called Self-Consistent RPA (SCRPA) equations |2(ll27ll28j |. extended 
to include higer configurations. In case the Xx 1 \ 2 \' 1 x' 2 amplitudes are dropped in eq. jJ53)l , eq.lO reduces to something 
similar to what has become known as renormalized RPA (r-RPA) p9l |. The main difference seems to come from the 
fact that here ea. 1(45(1 serves to determine the occupation matrix n° aa , whereas in r-RPA ea. 1(45(1 is used to establish 
the single particle basis. It should be interesting to investigate this relation more in detail in the future. 



V. SUMMARY 



Necessary conditions that the spurious state associated with the translational motion and its double-phonon state 
have zero excitation energy in extended ERPA (ERPA) were investigated using the small amplitude limit of the time- 
dependent density-matrix theory (STDDM). The reason why STDDM was used is that it has a quite general form of 



11 



the ERPA kind based on the HF ground state. In the case of the single spurious state it is found that ERPA which 
keeps all components of the one-body amplitudes gives the spurious state at zero excitation energy. This does not 
depend on approximations for the two-body amplitudes as long as they are properly antisymmetrized. For example, 
ERPA with only the 2p-2h and 2h-2p components of the two-body amplitudes preserves this property of the single 
spurious state. In the case of the double spurious state, all components of the two-body amplitudes are found necessary 
to yield the mode at zero excitation energy. Of course, no truncation in single-particle space should be made in both 
cases. The Kohn theorem for the single and double Kohn modes and the continuity equations for transition densities 
and currents were also investigated and found to hold under the same conditions as those necessary for the spurious 
states. It was pointed out that STDDM inherently has asymmetry and non-hermiticity, although it conserves various 
physical properties as mentioned above. A formulation of ERPA with hermiticity was also presented using TDDM, 
in which it was discussed that a three-body correlation matrix needs to be included in the description of ground-state 
correlations. The investigations in this work were performed for the spurious translational motion. It seems, however, 
clear that analogous considerations can be made for any spontaneously broken symmetry. An interesting case could 
be the coupling of quark-antiquark to the four quark sector using, e.g. the Nambu-Jona-Lasinio model [20]. In this 
case one knows that chiral symmetry is spontaneously broken |3(1 ] and therefore in the chiral limit a double Goldstone 
mode (two pions) should appear. In the case of finite current quark masses analogous equations to those yielding the 
Kohn modes considered here should exist, actually well known as the Gellmann-Oakes-Renner relation [3l|. 



APPENDIX A 



When ip a is chosen to be an eigenstate of the mean field hamiltonian fea. l|47(l ^. eas. l|45(l and 146(1 become 

(e«> - e a )n° aa , = £ (CS lWAa (aA 3 | W |A 1 A 2 ) - C^ XlXa (A 1 A 2 | W |a / A 3 )) (Al) 

Ai A2 A3 

( e a' + e /3' ~ e a — e /3)Ca/3a'/3' = B a/3a'0' + ^aPa'P' + Hapa'p'i (A2) 



where 



B % a >{]> = ( a i A 2MA 3 A 4 } a [(<5 qAi -n° aXi )(5 0X2 ~ n° 0X2 )nl 3a ,nl i0 , 

Ax A2 A3 A4 

- <xAMs«>-nl 3a ,)(5 Xi p,n Xi p r )}, (A3) 

PaP<x>0> = J! ( A 1 A 2|w|A 3 A 4 )[((5 q a 1 (5/3A2 - <5aAi"-°A 2 _ "■qA 1 '5/3A2)Ca3A4q'/3' 
Ai A2 A3 A4 

- (S\ 3a >6 Xil 3> -6 X3a ,n° Xi p, -n° X3a ,8 Xi0 >)Cl pXlX2 ], (A4) 

H a/3a>l3> = Y ^ ^ 2 ' V \ ^ 3 ^ A ^ aXl (™W C \4.i3\ 2 0' ~ n \ 3 /3' C \ i /3\ 2 a' + C i 3 A 4 /3a'A 2 /3' ) 



Ai A2 A3 A4 



+ ^A 2 ( n A 4 /3'CA3aAia' ~ n A 4 a' ^A 3 qA 2 /3' + C Xi \ 3 u0' \ 2 a' ) 
0a'A 3 l,"QAi O A4^A2/5' — n 0\i ° A 4 aA 2/ 3' + qA 4 /3AiA2/3V 

— Sp' Xi (n°p X2 C X3aXia , ~ n^ X2 C Xs p Xia , + Cp X3aX2Xia ,)]. (A5) 

The three-body correlation matrix C^g ,g, , is also included in ea. (|A5|) . The equation for C^ 7a ,g,^, is obtained by 
neglecting four-body amplitudes and becomes 

— J J* — TJ* -4- J J* 

u a' f3'~f' af3~f ^ a' f3 ln f' f3~fa T u a' 0' '7' a7/3 

— K«'/3'7'q/37 — Va'(3<yip~fa + K*'/3'7'q7/3j (A6) 



where 



£W'/3'7' = " [^ X Ma' P')A{n° Xl C% X2Y ~ n$ Xi C° 7X27 , + n° aXl C° g ^,) 

Ai A2 

+ (\1\2\v\a' d')C% lXiX2l ,l (A7) 



12 

Ai A2 A3 

+ c °f}\ 1 \ 2 c °\ 3 f3>y - Ca/3A 1 /3'C , 7A3A 27 ' + a11 other exchange terms). (A8) 
Equations for correlation matrices of higher ranks may be formulated according to the truncation rules given in 

APPENDIX B 



The matrix elements of eas.l|57 |l -(|6U |l are explicitly shown below. 

Si(a'a:M) = n° x ,J a , x ~ n° a , x S a x> , (Bl) 

Tx(a'a : AiA 2 AiA' 2 ) = C A , A , aA2 <5 a 'A 1 + C a '\' 1 \ 1 \* S aX l - C a , A , AlA2 5 a , A ' 1 - Cy 1 x' 2 o l x 1 s a'X 2 , (B2) 

T 2 (a'/3'a/3 : AA') = C^^a'A + CQ'pxafyx' - C°, , X/3 8 a y - C x , a , al3 8^x (B3) 



S 2 (a'/3'a/3 : A 1 A 2 A' 1 A' 2 ) = W^K^/j) + C°, A , Q/3 ) - A(6 aX[ 6 px ,)(A(n a , Xi n% X2 ) + C° a , fJ , XlX2 ) 

+ F{a'f3'a(3 : AiA 2 A' 1 A' 2 ) - F{a'(3'a(3 : AiA 2 A 2 A' 1 ) 

- F(a'P'pa : AiA 2 A' 1 A' 2 ) + F(a'f3'fja : X 1 X 2 X' 2 X' 1 ) 

- F(A' 1 A 2 AiA 2 : a(3a'p') + F(X' 1 X' 2 X 2 X 1 : apol 0) 

+ F(A' 1 A 2 AiA 2 : a/3/3' a') - J F 1 (A' 1 A' 2 A 2 Ai : app'a'), (B4) 

where 

F(a'f3'a(3 : AiA 2 A' 1 A' 2 ) = ^ [^(n°, Al n^ A Jn A ^ + n^C^,^ + n°, Al C°, A , Aa/3 - n£, Al C^ A , Aa/3 

_ n a'\ 2 ^0'\' 2 \ 1 p + n %X 2 < ^a'\' 2 \ 1 p + C ! a'l3'\' 2 \ 1 \ 2 l3^ ( B5 ) 

The three-body correlation matrix is also included in ea. (|B5|) . The matrix elements A,B,C, and D are given in the 
following. For simplicity, terms containing , g , , are not shown. They appear in B, C, and D: Wherever there is 

a term containing mP a a , C% a , , , there exists a corresponding term with C^g ,p, ,. The expressions for A,B,C, and 
D are not unique and their symmetry properties are not necessarily apparent. Equations l|45|l . (|46(l . and (|67J) allow 
us to take other expressions and guarantee symmetry properties. 

A(a'a : XX') = (e a > - e a )(8x a 'n°y a - 8y a n° a , x ) 

+ J2[( 7 8\v\a\)(A(n° x ,y a , s ) + C Q x , a , j5 ) + (X'a'\v\S^(A(n° x n Sa ) + C% aX ) 

- (X'8\v\a 7 ) A (n lX n a , s + C° a , lSX ) - ( ia '\v\8X) A {n xll n» a + C° x , Sia )} 

- ]T(<a'7MWA'aC° 7 , A7 + la8\v\ai)8x a 'Cl, lll& \ 

B(a'[3'aP : AA') = ^{(A' 7 | V |a/3) A [^«, A 4 7 ) + C ^,^) + {a' P' \v\X 7 ) A [A{n° x , a n° p ) + C° x , ja0 ]} 

7 

+ H(a'[3'af3 : XX') - H(f3'a'a/3 : XX') + H*(apa'p' : X'X) - H*{[3aa'(3' : A' A) 
+ I(a'f3'a[3 : XX') - I{a'{3'/3a : XX') + I*{aPa'P' : X'X) - I*(af3/3'a' : X'X), 
C{a'a : AiA 2 A' 1 A' 2 ) = i?(A' 1 A' 2 AiA 2 : aa'), 

where 

H(a'f3'al3 : XX') = -8 Xa <{(e a + ep - e a > - ep>)C x , B , a0 

+ Y(( 7 8\v\aP)(A(n x , y n a , s ) + C x ,^ s ) 

- [(lS\v\ ai ') A (n xll n% 5 n Y[3 + n° x ,^, YSp - n% y C° x , Y50 ) 

-yS-y' 



(B6) 



(B7) 
(B8) 
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+ (jS\v\a 7 ')n° yp C° x 

- ( r yS\v\^y')A(nx'ynl, s n° la + n% 7 C^ YSa ~ np, y Cx> -,>*<*) 
+ (P'jlvlS^Ainl^n^ + n ya C° sx ,^ - n%C? v<ry ) 

+ (P'7\v\S 7 ')n xri C%^}}, (B9) 

I(a'P'aP : AA') = £{[<7*l*MU(n|U4* n v,J + <' 7 C °'A^ " 

+ ( 7 S\v\aX)(n o x , C o a ,^ 6 + n^C^)] 

- (X'-y\v\aS) A [A(n alx n ^)n% + n° sx C a ,^ + n%C° alp , Xl 

The matrix D is given as 

D(a'f3'af3 : XiX 2 X[X' 2 ) = (e Q / + ep> - e a - ep)S 2 {d ft' a/3 : XiX^X^) 

+ {X[X' 2 \v\^)A(A(n° a , Xi n% X2 ) + C° a , p , XlX2 ) + (a' 0' ^X^) A (A(nl, a n° x , p ) + C° x , x , a0 ) 

+ J(a'f3'a/3 : X 1 X 2 X' 1 X' 2 ) + J*{a(3a'(3' : A' 1 A 2 AiA 2 ) 
+ K(a'{3'a{3 : X 1 X 2 X' 1 X , 2 ) + L{a'f3'a(3 : AiA^A!,) 

+ all other exchange terms of K and L, (BH) 

where 

J(a'(3'ap : A 1 A 2 A' 1 A 2 ) = -A(6 Xia ,6 X20 ,){J2[(^\v\a(3) - J2(W\v\a<y')n° 7 , 

7 5 7' 

- (7S\v\fr')n%))}(A(nl,y x , s ) + C X , X ,^ S ) 

][(lS\v\ay)(nl, p C% 2lS - n° x ,pC% il5 - n° x ,^, x , ps + n^C^) 



- hS\v\^)(n x[a C° YX ^ s - n° x , a C%^ s - n x ,^, x , aS + n x ,^ ( , x , aS )]}, (B12) 
K(a'f3'a{3 : A 1 A 2 A' 1 A 2 ) = ^'{^I^M^M^a^a, + n^C^J 

7<5 

+ W'\v\X 2 5) A {n%n\ a n Q Kl + n%C° x , x , ai ) 
+ all other exchange terms] 

- ^[(75|H«7')AK i7 4,n° A ^n°, A2 + n^ 7 n°, 5 C°, y/3A2 + C^C^jJ 

7^7' 

- ( 1 P'\vW6) A (nl, a n Q x ,y Yp n° SX2 + n^^C^ + ^C^) 

+ all other exchange terms]}, (B13) 

L{cipal3 : A!A 2 AiA 2 ) = - £[(a'7MAiA 2 > (nj^n^n^ + "^C^a/j) 

7 

+ all other exchange terms] 

+ £[(a'7MAi^(™V°A;XA 2 7^ + n %^ n % C KK-f + C l'S^pC\ x , ai ) 

+ all other exchange terms]. (B14) 

Finally we discuss a relation between ea. (|56[) and a set of the STDDM equations (eqs.lO and QJ). When the ground 
state is approximated by the HF one, eas. (jBl|) - (|B4|) become 

Si(a'a : XX') = (f a - f a <)6 aX ,5 a , x , (B15) 
S 2 (a'[3'af3: A 1 A 2 A' 1 A 2 ) = A(6 ayi 5 0K )A(6 a , Xl d , X2 )F°, , a0 , (B16) 
Ti(a'a : AiA 2 A' 1 A 2 ) = T 2 {a' (3' a(3 : XX') = 0, (B17) 
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where 

F a'/3'a0 = fafpfa'fff' ~ fafpfa'fp'- (B18) 

Equations (|BB|) - (|B8|) and I|B11(I become the following: 

A(a'a : AA') = {{e a , - e a )5 aX ,5 a , x + (\'a'\v\a\) A (f a , ~ f a )](f y - fx), (B19) 

B(a'[3'aP : AA') = -[(f a fpfp> + fcJpfpW P'\v\a/3} A 5 a , x - (f a f fj f a > + f a f f 3f a >)(\'a'\v\a/3) A 6 p , x 

- {fpfa'fp' + ff3f a >ff3>)(a'l3'\v\\[3) A 6 a x> 

+ (Uf«U> + fJ a >fpi){al3'\v\\a) A 5 y](jx, - fx), (B20) 

C{a'a : AiA 2 AiA' 2 ) = B(AiA£AiA 2 : aa'), (B21) 

D(a'P'a(3 : AiA 2 A' 1 A' 2 ) = F Xi x 2 x' 1 x' 2 {(^ a ' + e P' £p)A{8 a x' 1 S P x' 2 )A^x 1 8^x 2 ) 

+ (1 - fa' - f0')(a'P'\v\\i\2)AA(6 aK 6 x 2 ) 

- (1 - fa - f(j){X' 1 X 2 \v\a(3) A A{5 a >x 1 8(3>x 2 ) 

+ {fa - fa')[(a' \[\v\a\i) A 8/ 3 ,x 2 8px 2 + (a'X 2 \v\aX 2 ) A S l3 >x 1 S[3x' 1 

- (a' X' 1 \v\aX 2 ) A Sp' Xl Spx' 2 ~ {a' X' 2 \v\aXi) jtSp* x 2 Spx[] 

+ (fp - fpW'X' x \v\^i)A$ a '\ 2 5 a \ 2 + (P'X'M^2) A 5 a -x 1 5ay 1 

- {P' X' l \v\(3X 2 ) A 8 a >x 1 5 a x' 2 - {P'\' 2 \v\(3\i) A 5 a >\ 2 8 a \^] 

- {fa - fp')[{f3'X' 1 \v\a\\) A 8 a 'x 2 8p\' 2 + {f3'X' 2 \v\aX 2 ) A 8 a 'x 1 5i3x' 1 

- (P'X' 1 \v\aX 2 ) A S a >x 1 8px 2 ~ {(i'X' 2 \v\aX 1 ) A 8 al x 2 8 f 3X' 1 ] 

- {fa - / / 3')[(/3'A' 1 |u|aAi}A^a'A 2 ^ + {f3'X' 2 \v\aX 2 ) A 8 a 'x 1 8 l 3x' 1 

- {P'X'MaX^^x^px^iP'X'MaX^AS^xJpx^}- (B22) 

If S\x and S 2 X which appear in the equation for X, that is, Bx + DX = ujS 2 X, are replaced by x and X , respectively, 
the equation for X is equivalent to eq.QJ. However, the replacement S\x — > x and S 2 X — > X in Ax + CX = u>S\X 
cannot give eq.JSJ because of the symmetric coupling between x and X. Since the expression for C (ea. (IB8|) is not 
unique as mentioned above, we can always take an expression for C which leads to the same coupling matrix as in 
eq.© in the HF limit. Such an expression for C is the following: 

C(a'a : AiA 2 A' 1 A' 2 ) = F x 1 x 2 x< 1 x< 2 (( a 'A 2 MAiA^^a - ( a 'X[ \v\Xi X 2 ) A 8x>. 2 a 

- (KX'MaX^ASx.a' + {X' l X' 2 \v\aX l ) A 8x 2a -). (B23) 
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